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Abstract 
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1 Introduction 

Many problems in the optimization theory and other fields of mathematics ap- 
pear to be linear over semirings with idempotent addition (the so-called idempotent 
superposition principle [0 , which is a natural analog of the well-known superposition 
principle in Quantum Mechanics). The corresponding approach is developed system- 
atically as Idempotent Mathematics or Idempotent Analysis, a branch of mathematics 
which has been growing vigorously last time (see, e.g., IIT|-[§]). 

One of the most important examples of idempotent semirings is the set Mmax = 
M U {— 00} equipped with operations © = max, = + (see section |). In general, 
there exists a correspondence between interesting, useful, and important constructions 
and results concerning the field of real (or complex) numbers and similar constructions 
dealing with various idempotent semirings. This correspondence can be formulated in 
the spirit of the well known N. Bohr's correspondence principle in Quantum Mechan- 
ics; in fact, the two principles are intimately connected (see [|, ^, |^ and sections ^, 
H and ^ below). We discuss idempotent analogs of some basic ideas, constructions, 
and results in traditional calculus and functional analysis; also, we show that the 
correspondence principle is a powerful heuristic tool to apply unexpected analogies 
and ideas borrowed from different areas of mathematics (see, e.g., |]l|-|^). 

The theory is well advanced and includes, in particular, new integration theory, 
new linear algebra, spectral theory, and functional analysis. Its applications include 
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various optimization problems such as multicriteria decision making, optimization on 
graphs, discrete optimization with a large parameter (asymptotic problems), opti- 
mal design of computer systems and computer media, optimal organization of par- 
allel data processing, dynamic programming, applications to differential equations, 
numerical analysis, discrete event systems, computer science, discrete mathematics, 
mathematical logic, etc. (see, e.g., and references therein). 

In section |^ we give a short heuristic introduction into Idempotent Mathematics. 
Section ^ contains definitions of basic concepts of idempotent arithmetic and several 
important examples. In sections we consider the notion of linearity in Idempotent 
Analysis and indicate some of its applications to idempotent linear algebra. 

Due to imprecision of sources of input data in real-world problems, the data usu- 
ally come in a form of confidence intervals or other number sets rather than exact 
quantities. Interval Analysis (see, e.g., |T^-|T^) extends operations of traditional cal- 



culus from numbers to number intervals, thus allowing to process such imprecise data 
and control rounding errors in computations. To construct the analog of calculus of 
intervals in the context of optimization theory and Idempotent Analysis, we develop 
a set- valued extension of idempotent arithmetic (see section |]). 

The interval extension of an idempotent semiring is constructed in sections |^ 



and |I0| The idempotent interval arithmetic appears to be remarkably simpler than its 
traditional analog. For example, in the traditional interval arithmetic multiplication 
of intervals is not distributive with respect to interval addition, while idempotent 
interval arithmetic conserves distributivity. 

A simple application of interval arithmetic to idempotent linear algebra is dis- 
cussed in section |12|. We stress that in the traditional Interval Analysis the set of all 
square interval matrices of a given order does not form a semigroup with respect to 
matrix multiplication: this operation is not associative since distributivity is lost in 
traditional interval arithmetic. On the contrary, in the idempotent case associativity 
is conserved. 

Two properties that make the idempotent interval arithmetic so simple are mono- 
tonicity of arithmetic operations and positivity of all elements of an idempotent semir- 
ing. In general, idempotent interval analysis appears to be best suited for treating 
the problems with order-preserving transformations of imprecise data. We stress that 
this construction provides another example of heuristic power of the idempotent cor- 
respondence principle. 

2 Dequantization and idempotent 
correspondence principle 

Let M be the field of real numbers and M+ be the subset of all non-negative numbers. 
Consider the following change of variable: 

u w = h\nu, 
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where u G IR+ \ {0}, h > 0; thus u = e^^^, w G M. Denote by the additional element 
— oo and by S the extended real line M U {0}. The above change of variable has a 
natural extension to the whole S by -D/j(0) = 0; also, we denote Dh{l) = = 1. 

Denote by Sh the set S equipped with the two operations ®h (generalized addition) 
and Qh (generalized multiplication) such that Dh is a homomorphism of {M+,+,-} 
to {S, ©/i, Qh}- This means that Dh{ui + U2) = Dh{ui) ®h Dh{u2) and Dh{ui ■ U2) = 
Dhiui) Qh Dh{u2), i.e., wi Qh W2 = wi + W2 and wi Qh W2 = h\n{e'"^/'' + It is 

easy to prove that w\ Qh 102 maxjiwi, ^2} as h ^ 0. 

Denote by Mmax the set S* = M U {0} equipped with operations Q = max and 
= +, where = — 00, 1 = as above. Algebraic structures in and Sh are 
isomorphic; therefore Mmax is a result of a deformation of the structure in M_|_. 

We stress the obvious analogy with the quantization procedure, where h is the 
analog of the Planck constant. In these terms, M+ (or M) plays the part of a 'quantum 
object' while Mmax acts as a 'classical' or 'semi classical' object that arises as the result 
of a dequantization of this quantum object. 

Likewise, denote by Mmin the set M U {0} equipped with operations © = min 
and = +, where = +00 and 1 = 0. Clearly, the corresponding dequantization 
procedure is generated by the change of variables u w = —hlnu. 

Consider also the set M U {0, 1}, where = —00, 1 = +00, together with the 
operations © = max and = min. Obviously, it can be obtained as a result of a 
'second dequantization' of M or M+. 

The algebras presented in this section are the most important examples of idempo- 
tent semirings (for the general definition see section ^) , the central algebraic structure 
of Idempotent Analysis. The basic object of the traditional calculus is a function de- 
fined on some set X and taking its values in the field M (or C) ; its idempotent analog 
is a map X ^ S, where X is some set and S = Mmin, Mmax, or another idempotent 
semiring. Let us show that redefinition of basic constructions of traditional calculus 
in terms of Idempotent Mathematics can yield interesting and nontrivial results (see, 
e.g., [Q-IQ for details and generalizations). 

Example 2.1. Integration and measures. To define an idempotent analog 
of the Riemann integral, consider a Riemann sum for a function (p{x), x E X = [a,b], 
and substitute semiring operations © and for traditional addition and multiplication 
of real numbers in its expression (for the sake of being definite, consider the semiring 

Mmax) ■ 

N N 
J2vix^)■A, ^ 0(/?(Xi)0A, = max {(f{xi) + Aj), 

1=1 1=1 ~ 

where a = xq < Xi < ■ ■ ■ < x^ = b, Ai = Xi — i = 1, . . . , N. As maxj Aj 0, 
the integral sum tends to 

/ ^{x) dx = sup (p{x) 

J X x&X 
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for any function yj: X — > Mmax that is bounded. In general, the set function 

m^{B) = sn]np{x), B C X, 

is called an M.^aa^-'i^easure on X; since m^iXJ^Ba) = sup^m^{Ba), this measure is 
completely additive. An idempotent integral with respect to this measure is defined 
as 

/ ijj{x) dm^ = / ^/j{x) Q (p{x) dx = sup {ijj{x) + (p{x)). 

Jx Jx xf^X 

Example 2.2. Fourier-Legendre transform. Consider the topological 
group G = M". The usual Fourier-Laplace transform is defined as 

ip{x) ^ ^(0 = / e^«-^¥p(x) dx, 
Jg 

where exp(z^-a;) is a character of the group G, i.e., a solution of the following functional 
equation: 

f{x + y) = f{x)fiy). 
The idempotent analog of this equation is 

f(x + y) = fix)Qfiy) = fix) + f{y). 

Hence 'idempotent characters' of the group G are linear functions of the form x i-^ 
i ■ X = ^iXi + • • ■ + inXn- Thus the Fourier-Laplace transform turns into 

(p{x) I— > (p{^) = ^ ■ X Q (p{x) dx = sup ■ X + (p{x)). 
Jg xeG 

This equation differs from the well-known Legendre-Fenchel transform (see, e.g., [p^] ) 
in insignificant details. 

These examples suggest the following formulation of the idempotent correspon- 
dence principle [Q: 

There exists a heuristic correspondence between interesting, useful, and 
important constructions and results over the field of real (or complex) num- 
bers and similar constructions and results over idempotent semirings in the 
spirit of N. Bohr's correspondence principle in Quantum Mechanics. 

So Idempotent Mathematics can be treated as a 'classical shadow (or counterpart)' 
of the traditional Mathematics over fields. 

3 Idempotent semirings: Basic definitions 

Consider a set S equipped with two algebraic operations: addition © and multipli- 
cation 0. The triple {S, ©, 0} is an idempotent semiring if it satisfies the following 
conditions (here and below, the symbol denotes any of the two operations 0, 0): 
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• the addition © and the multiphcation are associative: x-k^y-kz) = {x-ky) -kz 
for all x,y,z E S; 

• the addition © is commutative: x Q) y = y Q) x for aA\ x, y & S; 

• the addition © is idempotent: x © a; = x for all a; G S'; 

• the multiphcation is distributive with respect to the addition ©: xQ{yQ)z) = 
{x &y) (B {x Q z) and {x ® y) Q z = {x Q z) ® {y Q z) for all x,y,z e S. 

In the rest of this paper we shall sometimes drop the word 'idempotent' when the 
context is clear. 

A unity of an idempotent semiring S is an element 1 G S* such that for all s G 5* 

l(Dx = xQl = x. 

A zero of an idempotent semiring S is an element G S* such that 7^ 1 and for 
all a; G S* 

0©a; = a;, O0x = a;0O = O. 

It is readily seen that if an idempotent semiring S contains a unity (a zero), then 
this unity (zero) is determined uniquely. 

A semiring S is said to be commutative if x Q y = y Q x for all x,y & S. A 
commutative semiring is called a semifield if every nonzero element of this semiring 
is invertible. It is clear that M^ax and Mmin are semifields. 

Note that different versions of this axiomatics are used, see, e.g., and some 

literature indicated in these books and papers. 

The addition © defines on an idempotent semiring 5* a partial order: x ^ y iS 
X ® y = y. We use the notation x -< y ii x ^ y and x ^ y. li S contains a zero 0, 
then is its least element with respect to the order The operations © and are 
consistent with the order ^ in the following sense: ii x ^ y, then x -k z ^ y -k z and 
z -k X ^ z -k y for all x, y, z ^ S. 

An idempotent semiring S is said to be a- complete if for any subset {xa} C S, 
including 0, there exists a sum 0{xq,} = 0^ Xa such that (0^ Xa)Qy = ^a^x^Qy) 
and yQi^^Xa) = 0Q,(l/0a;a) for any y & S. An idempotent semiring S containing a 
zero is said to be b-complete if the conditions of a-completeness are satisfied for any 
nonempty subset {xa} C S that is bounded from above. Any 6-complete semiring 
either is a-complete or becomes a-complete if the greatest element 00 = sup 5* is 
added; see P, ^ for details. 

Note that 0^ Xq, = sup{xq,}; in particular, an a-complete idempotent semiring 
always contains the zero = 00. 

An idempotent semiring S does not contain zero divisors if x ?/ = implies 
that X = or y = for all X, ?/ G S. An idempotent semiring S is said to satisfy 
the cancellation condition if for all x, y, z & S such that x 7^ it follows from 
X Q y = X Q z OT y Q X = z Q x that y = z. Any idempotent semiring satisfying 
the cancellation condition does not contain zero divisors. A commutative idempotent 
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semiring S is said to be idempotent semifield if every nonzero element of 5* is invertible; 
in this case the cancellation condition is fulfilled. 

An idempotent semiring S is said to be algebraically closed if the equation x" = y, 
where a:" = x Q x d) ■ ■ ■ Q x {n times), has a solution for all y G S" and n G N; 
an idempotent semiring S with a unity 1 satisfies the stabilization condition if the 
sequence x" © y stabilizes whenever x ^ 1 and y |jl5|, |16[. Note that in |]l6l the 



property of algebraic closedness was incorrectly called 'algebraic completeness' due to 
translator's mistake. 

The most important examples of idempotent semirings are those considered in 
section 0. We see that M^ax is a 6-complete algebraically closed idempotent semi- 
field satisfying stabilization condition. The idempotent semiring Mmin is isomorphic 
to Mmax- Note that both Mmax and Mmin are linearly ordered with respect to the cor- 
responding addition operations; the order =^ in Mmax coincides with the usual linear 
order ^ and is opposite to the order ^ in Mmin- 

Consider the set Mmax = MmaxUjoo} with operations ©, extended by oo©x = oo 
for all X G Mmax, oo0x = ooifx7^O and oo = 0. It is easily shown that this set 
is an a-complete idempotent semiring and oo is its greatest element. 

Let {5*1, ^2, . . .} be a collection if idempotent semirings. There are several ways 
to construct a new idempotent semiring derived from the semirings of this collection. 

Example 3.1. Suppose 5* is an idempotent semiring and X is an arbitrary set. 
The set A1(X; S) of all functions X — S* is an idempotent semiring with respect to 
the following operations: 

{f®9){x) = f{x)®g{x), {fQg){x) = f{x)Qg{x), x e X. 

If S contains a zero and/or a unity 1, then the functions o{x) = for all x E X, 
e{x) = 1 for all X G X are the zero and the unity of the idempotent semiring A4{X; S). 
It is also possible to consider various subsemirings of M{X; S). 

Example 3.2. Let Si be idempotent semirings with operations ©j, 0j and zeros 
Oj, i = 1, . . . ,n. The set 5* = (5*1 \ {Oi}) x • ■ ■ x (Sn \ {On}) U is an idempotent 
semiring with respect to the following operations: x-ky = (xi, . . . , x„) * . . . , = 
yn); the element is the zero of this semiring. 

Note that the set S* = S'l x ■ ■ ■ x is also an idempotent semiring with respect 
to the same operations; its zero is the element (Oi, . . . , 0„). 

Notice that even if primitive semirings in examples and |3.2| are linearly ordered 
sets with respect to the orders induced by the correspondent addition operations, the 
derived semirings are only partially ordered. 

Example 3.3. Let S be an idempotent semiring and MatmniS) be a set of all 
5*- valued matrices. Define the sum © of matrices A = (aij), B = (btj) G Matmn('S') as 
AQ) B = {aij © bij) G Matmn(S'), and let ^ be the corresponding order on Matrrtn(5'). 
The product of two matrices A G Matim{S) and B G Matmn('S') is a matrix AB G 
Matin{S) such that AB = Oifc bkj). The set Matn„(S') of square matrices is 

an idempotent semiring with respect to these operations. If is a zero of S, then the 
matrix O = (ojj), where Oij = 0, is the zero of Mat„„(S'); if 1 is a unity of S, then 
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the matrix E = (Sij), where 6ij = 1 if i = j and 6ij = otherwise, is the unity of 
Mat„„(5). 

Many additional examples can be found, e.g., in 

4 Idempotency and linearity 

Now we discuss an idempotent analog of a linear space. A set V is called a 
semimodule over an idempotent semiring S (or an S'-semimodule) if it is equipped with 
an idempotent commutative associative addition operation ©y and a multiplication 
0y: S X V ^ V satisfying the following conditions: for all X, fi & S, v, w E V 

• (A /i) 0y = A 0y (/i 0y v)] 

• XQv {v ®v w) = (A 0y v) ®v (A Qv w); 

• (A © /i) 0y = (A Qv V) ©y (/i 0y v). 

An S'-semimodule V is called a semimodule with zero if 5* is a semiring with a zero 
G S* and there exists a zero element Oy & V such that for all v E V, \ E S 

• Oy ©y t; = f ; 

• A 0y Oy = 0y W = Oy. 

Example 4.1. Finitely generated free semimodule. The simplest S'- 
semimodule is the direct product S" = { (ai, . . . , a„) \ aj E S, j = 1, . . . ,n}. The set 
of all endomorphisms 5"" —>■ coincides with the semiring Mat„„(5') of all 5- valued 
matrices (see example |373|) . 

Example 4.2. Matrix semimodule. Take some c e S, A e MatmniS). The 
product c Q A is defined as the matrix (c aij) G Matmn{S). Then the set of all 
S"- valued matrices of a given order Matmn(5') forms a semimodule under addition © 
and multiplication by elements of S. 

Example 4.3. Function spaces. An idempotent function space J^{X;S) is a 
subset of the set of all maps X ^ S such that if f{x), g{x) G J-'{X; S) and c E S, then 
(/ © 9)ix) = fix) © g{x) G J^(X; S) and (c © /)(x) = c /(x) G J^(X; S); thus an 
idempotent function space is another example of an S'-semimodule. If the semiring 5' 
contains a zero element and J-'{X; S) contains the zero constant function o{x) = 0, 
then the function space J-'{X] S) has the structure of a semimodule with zero o{x) 
over the semiring S. 

Recall that the idempotent addition defines a partial order in semiring S. An 
important example of an idempotent functional space is the space B{X; S) of all 
functions X ^ S bounded from above with respect to the partial order =<; in S. 
There are many interesting spaces of this type including C{X; S) (a space of continuous 
functions defined on a topological space X), analogs of the Sobolev spaces, etc. (see, 
e.g., 0-0 for details). 
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According to the correspondence principle, many important concepts, ideas and 
results can be converted from usual functional analysis to Idempotent Analysis. For 
example, an idempotent scalar product in B{X] S) can be defined by the formula 

{ip,ip) = / ip{x) Q ipi^x) dx, 
Jx 

where the integral is defined as the sup operation (see example Notice, however, 
that in the general case the ordering ^ in is not linear. 

Example 4.4. Integral operators. It is natural to construct idempotent 
analogs of integral operators of the form 

K : ip{y) {Kip){x) = J K{x, y) Lp{y) dy, 

where (p{y) is an element of a functional space J-'i{Y; S), {Kip){x) belongs to a space 
^2(^5 S) and K{x,y) is a function X x Y S. Such operators are homomorphisms 
of the corresponding functional semimodules. If 5 = Mmax, then this definition turns 
into the formula 

{K(p){x) = sup {K{x,y) + (fiy)). 

y& 

Formulas of this type are standard for optimization problems (see, e.g., jlTll). 



5 Idempotent superposition principle 

In Quantum Mechanics the superposition principle means that the Schrodinger 
equation (which is basic for the theory) is linear. Similarly in Idempotent Mathemat- 
ics the (idempotent) superposition principle means that some important and basic 
problems and equations (e.g., optimization problems, the Bellman equation and its 
versions and generalizations, the Hamilton- Jacobi equation) nonlinear in the usual 
sense can be treated as linear over appropriate idempotent semirings, see |T]1^0- 

The linearity of the Hamilton- Jacobi equation over Mmin (and Mmax) can be de- 
duced from the usual linearity (over C) of the corresponding Schrodinger equation 
by means of the dequantization procedure described above (in section |^). In this 
case the parameter h of this dequantization coincides with ih , where h is the Planck 
constant; so in this case h must take imaginary values (because h > see 0, ^ for 
details). Of course, this is closely related to variational principles of mechanics; in 
particular, the Feynman path integral representation of solution to the Shrodinger 
equation corresponds to the Lax-Olemik formula for solution of the Hamilton- Jacobi 
equation (for the latter see, e.g., |T8[] ). 

The situation is similar for the differential Bellman equation, see 0, |^ . 

B.A. Carre used the idempotent linear algebra to show that different optimiza- 



tion problems for finite graphs can be formulated in a unified manner and reduced to 
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solving Bellman equations, i.e., systems of linear algebraic equations over idempotent 
semirings. 

Discrete Bellman equation. It is well-known that discrete versions of the 
Bellman equation can be treated as linear over appropriate idempotent semirings. The 
following equation (the discrete stationary Bellman equation) plays an important role 
in both discrete optimization theory and idempotent matrix theory: 

X = AX®B, (1) 

where A G Matnn(5'), X,B ^ Matri;(S'); matrices A, B are given and X is unknown. 
Equation (|ID is a natural counterpart of the usual linear system AX = B m traditional 
linear algebra. 

Note that if the closure matrix A* exists, then the matrix X = A*B satisfies (|l]) 
because A* = AA* © E. It can be easily checked that this special solution is the 
minimal element of the set of all solutions to (|I]). 

Actually, the theory of the discrete stationary Bellman equation can be devel- 
oped using the identity A* = AA* © as an additional axiom (the so-called closed 
semirings; see, e.g., ). 



B. A. Carre []T^ also generalized to the idempotent case the principal algorithms of 
computational linear algebra and showed that some of these coincide with algorithms 
independently developed for solution of optimization problems; for example. Bellman's 
method of solving shortest path problems corresponds to a version of Jacobi's method 
for solving systems of linear equations, whereas Ford's algorithm corresponds to a 
version of Gauss-Seidel's method. 

We stress that these well-known results can be interpreted as a manifestation of 
the idempotent superposition principle. 



6 Idempotent matrix theory: Some results 

Matrix algebra and graph theory. Any square matrix A = (aij) e Mat„ri(5') 
specifies a weighted directed graph with n nodes such that aij G S' is the weight of 
the arc connecting ith node to jth one. On the other hand, let G be a directed graph 
with at most one arc between any two nodes, every arc of which is characterized by 
its weight that belongs to S. Then G can be described by a matrix A G Mat„,i(S'); in 
particular, nonexistent arcs are counted with weight 0, and nonzero diagonal entries 
an of matrix A correspond to loops. 

For any A G Mat„„(5) define A^ = E, A'' = AA''-\ k > 1. Let a'^^ be {i,j)th 
element of the matrix A'^; then it is easy to check that 

io=i, l^ii,...,ifc_i^n,ifc=j 

Consider a path of length k in the graph G, i.e. a sequence of nodes iq, . . . ,ik connected 
by k arcs with weights ajgj^, . . . , ai^_^i^. The weight of the whole path is defined to be 
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the product aj^jj ■ ■ ■ ai^_^i^- Thus af^ is the supremum of weights of all paths of 
length k connecting node Iq = i to node = j- 

Algebraic path problem. This well-known problem is formulated as follows: 
for each pair calculate supremum of weights of all paths (of arbitrary length) 
connecting node i to node j. If the semiring 5* under investigation is Mmin and arc 
weights are lengths in some metric, then to solve this problem means to find all shortest 
paths. Is 5* is {0, 1} and the corresponding directed graph depicts some relation R in 
the set {1, . . . ,n}, then to solve this problem means to find the transitive closure of 
R. 

It is evident that in terms of matrix theory the algebraic path problem is reduced 
to finding the matrix satisfying the formal expansion 



A* = E(BAQ)A^®--- = ^A''. 



k=0 

The matrix A* is called the closure of the matrix A. If the idempotent semiring 
S is not a-complete, this problem can be nontrivial since we cannot simply take the 
infinite sum. Below we shall discuss a sufficient condition for the existence of a closure, 
following work of B. A. Carre ||19|| . 

The matrix A = (aij) G Mat„„(S') is said to be definite [semi- definite) if 

C^joU ■ ■ ■ © ^ik-iik ^ ('^io*l ■ ■ ■ ^ijc-iik ^ 

for any path {ii, . . . ,ii) such that Iq = i^ (i.e., for any closed path). Obviously, every 



definite matrix is semi-definite. This definition is similar to that of B. A. Carre |jT9 
with the only difference: Carre considers an ordering that is opposite to ^. 

Theorem 1 (Carre) Let A he a semi- definite matrix. Then 

k n— 1 

®A' = ®A' 

1=0 1=0 

for k ^ n — 1, so the closure matrix A* = A^ exists and is equal to ©^Zq 

For the proof see Theorem 4.1]. The basic idea of the proof is evident: in the 
graph of a semi-definite matrix it is impossible to construct a path of arbitrarily large 
weight since the weight of any closed part of a path cannot be greater than 1. Thus 
there exists a universal bound on path weights, which makes possible truncation of 
the infinite series for the closure matrix. 

Spectral theory. The spectral theory of matrices whose elements lie in an 
idempotent semiring is similar to the well-known Perron-Frobenius theory of nonneg- 
ative matrices (see, e.g., [0, ||, |15|, p!6|). 

Recall that the matrix A = (aij) G Mat„„(S') is said to be irreducible in the sense 
of if for any 1 ^ i, j ^ n either a^j 7^ or there exist 1 ^ ii, . . . ,ik ^ n such that 
au-^ ■ ■ -Qai^j 7^ 0. In |jl5|, [1^ matrices with this property are called indecomposable. 



We borrow the following important result from |T^ (see also 
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Theorem 2 (Dudnikov, Samborskii) If a commutative idempotent semiring S with 
a zero and a unity 1 is algebraically closed and satisfies cancellation and stabiliza- 
tion conditions, then for any matrix A G Mat„„(5') there exist a nonzero 'eigenvector' 
V e Matni(5') and an 'eigenvalue' \ E S such that AV = XqV. If the matrix A is 
irreducible, then the 'eigenvalue' \ is determined uniquely. 

For the proof see [|TB], Theorem 6.2]. 



An apphcation of this result will be given in section 12 



Similar results hold for semimodules of bounded or continuous functions 



7 Correspondence principle for computations 

Of course, the (idempotent) correspondence principle is valid for algorithms as well 
as for their software and hardware implementations 0, 0. Thus: 

// we have an important and interesting numerical algorithm, then there 
is a good chance that its semiring analogs are important and interesting 
as well. 

In particular, according to the superposition principle, analogs of linear alge- 
bra algorithms are especially important. Note that numerical algorithms for stan- 
dard infinite-dimensional linear problems over idempotent semirings (i.e., for prob- 
lems related to idempotent integration, integral operators and transformations, the 
Hamilton- Jacobi and generalized Bellman equations) deal with the corresponding 
finite-dimensional (or finite) "linear approximations". Nonlinear algorithms often 
can be approximated by linear ones. Thus idempotent linear algebra is the basis of 
the idempotent numerical analysis. 

Moreover, it is well-known that algorithms of linear algebra are convenient for 
parallel computations; their idempotent analogs admit parallelization as well. Thus we 
obtain a systematic way of applying parallel computation to optimization problems. 

Basic algorithms of linear algebra (such as inner product of two vectors, matrix 
addition and multiplication, etc.) often do not depend on concrete semirings, as well 
as on the nature of domains containing the elements of vectors and matrices. Thus 
it seems reasonable to develop universal algorithms that can deal equally well with 
initial data of different domains sharing the same basic structure |H, H]; an example 



of such universal Gauss- Jordan elimination algorithm is found in ^0 . 

Numerical algorithms are combinations of basic operations with 'numbers', which 
are elements of some numerical domains (e.g., real numbers, integers, etc.). But every 
computer uses some finite models or finite representations of these domains. Discrep- 
ancies between 'ideal' numbers and their 'real' representations lead to calculation 
errors. This is another reason to deal with universal algorithms that allow to choose 
a concrete semiring and take into account the effects of its concrete finite representa- 
tion in a systematic way; see || for details and applications of the correspondence 
principle to hardware and software design. 



11 



8 Set-valued idempotent arithmetic 



Suppose S is an idempotent semiring and 5 is a system of its subsets. We shall 
denote the elements of 5 by x, y, . . . and define x*y = {t & S \ t = x-ky, x G x, y G y }. 
We shall suppose that S satisfies the following two conditions: 

• if X, y G and * is an algebraic operation in S, then there exists z G 5 such 
that z D X y; 

• if {xq,} is a subset of S such that n«^a 7^ then there exists the infimum of 
{xq,} in S with respect to the ordering C, i.e., the set y G 5 such that y C f]^ Xq, 
and z C y for any z E S such that z C f]^ Xq,. 

Define algebraic operations ©, © in 5 as follows: if x, y G S, then x^y is the 
infimum of the set of all elements z G 5 such that z D x^y. 

Proposition 1 The following assertions are true: 

1. S is closed with respect to operations ®, 0. 

2. The element x^y zs optimal in the following sense: suppose the exact values of 
input variables x and y lie in sets x and y, respectively; then the result of an 
algebraic operation x ^ y contains the quantity x-ky and is the least subset of S 
in S with this property. 

3. If the systems contains all one-element subsets of S, then the semiring {5*, ©, 0} 
is isomorphic to a subset of the algebra {S, 0, 0}. 

The proof is straightforward. 

In general, not much can be said about the algebra {S, 0, 0}, as the following 
example shows. 

Example 8.1. Let S = 2^ and x^y = x^y. In general, the set S with these 
'naive' operations 0, satisfies the above assumptions but is not an idempotent 
semiring. Indeed, let S be the semiring (Mmax \ {0})^ U {(0,0)} with operations 0, 
defined as in example [X^ . Consider a set x = {(0, 1), (1, 0)} G 5; we see that 
X ©X = {(0, 1), (1, 0), (M)}/ X and if y = {(1, 0)}, z = {(0, 1)}, then x0 (y © z) = 
{(1^), (2, 1)} ^ (x0y)©(x0z) = {(1, 1), (1, 2), (2, 1), (2, 2)}. Thus S with operations 
0, does not satisfy axioms of idempotency and distributivity. 

It follows that S should satisfy some additional conditions in order to have the 
structure of an idempotent semiring. In the next sections we consider the case when 
5 is a set of all closed intervals; this case is of particular importance since it represents 
an idempotent analog of the traditional Interval Analysis. 
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9 Weak interval extension of an idempotent 
semiring 

Let S be an idempotent semiring. A (closed) interval in S" is a set of the form 
x = [x, x] = {te5'|x^^^x}, where x, x G 5* (x ^ x) are said to be lower and 
upper bounds of the interval x, respectively. 

Note that if x and y are intervals in S, then xCyiffy^x=>>x^y. In 
particular, x = y iff x = y and x = y . 

Remark 9.1. Let x, y be intervals in S. In general, the set x T<r y is not an 
interval in S. Indeed, consider a set S* = {0, a, b, c, d} and let © be defined by the 
following order relation: is the least element, d is the greatest element, and a, b, and 
c are noncomparable with each other. If © is a zero multiplication, i.e., if x Q y = 
for all X, y & S, then S is an idempotent semiring without unity. Let x = [0, a] and 
y = [0, b]; thus x©y = {0, a, b, d}. This set is not an interval since it does not contain 
c although =4 c ^ d. 

Let S be an idempotent semiring. A weak interval extension I{S) of the semiring 
S is the set of all intervals in S equipped with the operations ©, ©, where for any two 
intervals x, y G I{S) x*y is defined as the smallest interval containing the set x^y. 

Proposition 2 I{S) is closed with respect to the operations ©, © and has the struc- 
ture of an idempotent semiring. Moreover, x^y = [x*y,X7iry] for a//x, y G I{S). 



Proof. Take t G x t^t y and let x G x, y G y be such that t = x-ky. By definition of 
interval, it follows that x ^ x ^ x and y ^ y ^ y- Since the operation -k is consistent 
with the order we see that x^y^x^iry^x^y. In particular, x*y ^ x^^ry, i.e., 
the interval [x^y,x*y] is well defined. It follows that x*y C [x*y,x^y]. 

Now let an interval z in 5* be such that x^y C z. We have x^y^x^yCz and 
x^yGx^yCz; hence z ^ x*y and x^y ^ z. Since x^y^x^yby the above, it 
follows that [x*y,XT<ry] C z. 

Thus the set I{S) with the operations ©, © can be identified with a subset of an 
idempotent semiring S x S (see example |3.2| ). Since x*y =^ x^^ry whenever x =^ x 
and y =^ y , I{S) is closed with respect to ©, ©; hence it is an idempotent semiring 
(a subsemiring oi S x S). □ 

Remark 9.2. Note that I{S) satisfies the third condition of proposition |l| only 
if the semiring 5* is 6-complete; nevertheless, the operations ©, © are well defined in 
the general case. 

Proposition 3 If S is an a-complete (b-complete) idempotent semiring, then I{S) is 
an a-complete (b-complete) idempotent semiring. 

Proof. Let 5* be an a-complete (6-complete) idempotent semiring and {x^} be 
a nonempty subset of I{S) (a nonempty subset of I{S) such that the set {xq} is 
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bounded in S). We claim that the interval [0^Xa) ©aXa] contains the set = { t G 
S I (Va)(3xQ, G Xq,) t = Xa } and is contained in every other interval containing iV; 
hence 0^^^" = [®a'^i®a^a\- Indeed, in the case of an a-complete semiring this 
statement is proved similarly to proposition 0. If S is 6-complete and the set {xa} is 
bounded from above, then the set {Xq,} is also bounded from above, i.e., there exists 
y E S such that x^ ^ Xq =^ ?/ for all a. Thus there exist 0q,Xq, G S and 0q,Xq, G S. 
Now the obvious adaptation of the above proof completes the argument. 

If S is a-complete and X C I{S) is empty, then by above definitions 0X = [0, 0] 

and for all y G I{S) y = (©^) ©y = [0,0]. A direct calculation shows 

that if X = {xq,} C I{S) is nonempty, then 

for any y G I{S). Thus /(S) is a-complete (6-complete) if S is a-complete (6- 
complete) . □ 
The following two propositions are straightforward consequences of proposition 

Proposition 4 If an idempotent semiring S is commutative, so is I{S). 



Proposition 5 // an idempotent semiring S contains a zero {unity 1), then the 
interval [0,0] ([1,1]) is the zero (unity) ofI{S). 



Proposition 6 // an idempotent semiring S has a zero and does not contain zero 
divisors, I{S) has no zero divisors as well. 

Proof. Let x, y G I{S) and x 7^ [0, 0], y 7^ [0, 0]. Because x ^ x, y ^ y , this means 
that X 7^ 0, y 7^ 0. Thus if z = x y, then z = x y 7^ 0, since there are no zero 
divisors in S. It follows that z 7^ [0,0]. □ 

Proposition 7 If S is algebraically closed and for any x,y & S , n & N the equality 
(x © y)" = x" © ?/" holds, then I{S) is algebraically closed. 

Proof. Suppose x" = x x ■ ■ ■ x = y. By proposition ^ we see that x" = y 
and x*^ = y. Let z G S' and z G S' be the solutions of these two equations. We claim 
that z and z can be chosen such that z ^ z, i.e., the interval [z, z] is well defined. 

Take z' = z©z; hence z ^ z'. Since in S z'" = (z©z)" = z"©z", z"" = y ©y = y. 
We see that [z,z']'' = [y,y] = y. ~ □ 

Remark 9.3. Note that the equality (x©?/)" = x"©?/" holds in any commutative 
idempotent semiring 5* satisfying cancellation condition (see, e.g., [T^|, assertion 2.1). 
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10 A stronger notion of interval extension 



We stress that in general a weak interval extension I{S) of an idempotent semiring S 
that satisfies cancellation and stabilization conditions does not inherit the latter two 
properties. Indeed, let x © z = y z, where z = [0, z] and -2 7^ 0; then z is a nonzero 
element but this does not imply that x = y since x and y may not equal each other. 
Further, let y = [0, y] 7^ [0, 0]; then the lower bound of x"0y may not stabilize when 
n — >• cx). 

Therefore we define a stronger notion of interval extension of an idempotent semir- 
ing S with a zero to be the set I(S') = {[x,y] \ x,y & S, ^ x ^ y, } U {[0,0]} 
equipped with operations ©, © defined as above. 

Note that this object may not be well-defined if the semiring 5* has zero divisors. 
Indeed, let x = [a;i,X2] G I{S) and y = [yi,y2] G i{S) be such that -< xi -< X2, 
-< 1/1 -< ?/2, 2^1 2/1 = 0, and X2 © 2/2 7^ 0; then x © y = [0, 0:2 © 2/2] ^^S). 

Throughout this section, we will suppose the interval extension 1(5) of an idem- 
potent semiring S to be closed with respect to the operations © and ©. To achieve 
this, it is sufficient to require that the semiring S contains no zero divisors. 

Theorem 3 The set I{S) is an idempotent semiring with respect to the operations © 
and © with the zero [0, 0] and does not contain zero divisors. It inherits some special 
properties of the semiring S: 

1. If S is a-complete (h-complete), then I(S') is a-complete (h-complete). 

2. If S is commutative, so is 1(5'). 

3. If 1 is a unity of S, [1,1] is the unity o/I(5). 

4- If S is algebraically closed and for any x,y E S , n E N the equality (x © y)"^ = 
x" © holds, then 1(5) is algebraically closed. 

5. If S satisfies the cancellation condition, so does 1(5). 

6. If S satisfies the stabilization condition, so does 1(5). 

Proof. Using proposition ||, it is easy to check that 1(5) is an idempotent 
semiring with respect to the operations ©, ©. This semiring has the zero element 
[0, 0] but contains no zero divisors by propositions ^ and ^ Propositions and |^ 
imply the first four statements. 

Suppose 5 satisfies the cancellation condition, x, y, z G 1(5), and z is nonzero. 
If X © z = y © z, then x © z = y © z and x © z = y © z; since z 7^ [0, 0] in 1(5), 
z 7^ and z 7^ 0, and it follows from the assumptions that x = [x, x] = [y,y] = y. If 
z © X = z © y, then x = y similarly. 

Suppose further that 5 satisfies the stabilization condition; by definition, y 7^ 
and y 7^ for any nonzero y G 1(5). Consider the sequence x" © y; stabihzation 
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holds in S for both bounds of the involved intervals and hence, by proposition ^, for 
the whole intervals as elements of 1(5'). □ 
Suppose S is an idempotent semiring; then the map l: S ^ I{S) defined by t{x) = 
[x,x] is an isomorphic imbedding of S into its weak interval extension I{S). If the 
semiring S has a zero but no zero divisors, then the map l: S ^ I(S') C I{S) is an 
isomorphic imbedding of S into its interval extension. In the sequel, we will identify 
the semiring 5* with subsemirings l{S) C I{S) or l{S) C I{S) C I{S) and denote the 
operations in I{S) or 1(5') by ©, 0. If the semiring S contains a unity 1, then we 
denote the unit element [1, 1] of I{S) or I{S) by 1; similarly, we denote [0,0] by 0. 



11 Cancellation and semifields 

We stress that in idempotent interval mathematics most of algebraic properties of an 
idempotent semiring are conserved in its interval extension. On the contrary, if S 
is an idempotent semifield, then the set 1(5) is not a semifield but only a semiring 
satisfying the cancellation condition. 

Recall that any commutative idempotent semiring 5 with a zero can be iso- 
morphically embedded into an idempotent semifield 5" provided that 5" satisfies the 



cancellation condition (see, e.g., |jT5|). If 5" coincides with its subsemifield gener- 
ated by 5, then 5* is called a semifield of fractions corresponding to the semiring 5". 
This semifield can be constructed as the quotient 5 x (5 \ {0})/ ~, where for any 
{x,y),{z,t)eSx{S\{0}) 

{x,y) ^ {z,t) <^ xQt = yQz, 

equipped with operations 

(x, y) © {z, t) = {{xQt)®{yQz),yQt), {x, y) (z, t) = {x Q z,y Q t). 

It is easy to see that these operations are defined in such a way that pairs {x,y) are 
treated as 'fractions' with x as 'numerator' and y as (nonzero) 'denominator'. This 
semifield has the zero element { (0, ?/) \ y ^ Q} and the unit element { (y, y) | y 7^ }; 
for every 'fraction' (x, y) representing a nonzero element of 5* its inverse element is 
given by the fraction (y, x). 

In the context of the traditional Interval Analysis a similar extension of the algebra 
of numerical intervals leads to the so-called Kaucher interval arithmetic p2i . In 
addition to usual intervals [x, y], where x ^ it includes quasiintervals [x,y\ with 
y ^ X, which arise as inverse elements for the former. 

Here we describe an idempotent version of Kaucher arithmetic. The following 
statement shows that in this case the semifield of fractions of interval extension 1(5") 
corresponding to an idempotent semiring 5* with cancellation condition has very sim- 
ple structure: it is isomorphic to the idempotent semifield {S \ {0}Y U {0}. 

Proposition 8 Suppose S is a commutative idempotent semiring with a zero 0, S 
satisfies the cancellation condition, and S is its semifield of fractions; then a semifield 
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of fractions corresponding to the interval extension I{S) is isomorphic to the semifield 
(5\ {0})2 U {0} (see example W^). 



Proof. It follows from theorem |^ that I(S') is a commutative idempotent semiring 
with the zero element = [0, 0] and satisfies the cancellation condition. Thus 1(5*) 
can be isomorphically embedded into its semifield of fractions. 

Define the map 0:1(5) x (1(5) \ {0}) ^ (5\ {0})^ U {0} by the rule (/'((x,y)) = 
(x y~^, X y^"*^). This map is surjective. Indeed, (0,0) = 0((O,y)) for any y 7^ 0; 
let us check that if a, 6 G S", a 7^ 0, 6 7^ 0, then there exist x, y G 1(5), y 7^ 0, such 
that (a, b) = 0((x, y)). We see that there exist ai, 02, &2 ^ S such that a = ai0a2 ^, 
b = biQb2^. Define x = ai 61 62, x = (ai 61 62) © (02 6?), y = 02 &i &2, 
y = (ai 62) © (0,2 Q bi Q 62); thus x =^ x, y ^ y and x y"^ = Oi 02 ^ = a, 
x0y-i = 61 062^ = 6. 

Since x = z Q iff x t = y Q z ioi any x,?/, z,t G 5, we see that 
(/)((x, y)) = 0((z,t)) iff (x, y) and (z,t) define the same element of the semifield of 
fractions corresponding to 1(5). Also, 

0(((x t) © (y z), y t)) = ((x y-i)) © (z t" '), (x y~')) © (z t"')) 

= 0((x,y))©0((z,t)), 
</.((x z, y t)) = ((x y-i) (z t" '), (x y"') (z t ')) 
= 0((x,y))0 0((z,t)). 

Thus the semifield of fractions corresponding to 1(5) is isomorphic to the idempotent 
semifield (5 \{0})2u{0}. □ 
The commutativity condition in this proposition is a natural one. Indeed, it can 
be proved that if each nonzero element of a 6-complete idempotent semigroup 5 has 
a multiplicative inverse, then 5 is commutative (see, e.g., f^). 

12 Application to linear algebra 

Suppose 5 is an idempotent semiring with a zero and a unity 1 and /(5) is its 
weak interval extension; then Mat„„(/(5)) is an idempotent semiring. If the interval 
extension 1(5) of the semiring 5 is well defined, then the same is true for Mat„„(I(5)). 
We shall denote the (common) unit element of these semirings by E. 

If A = (ajj) G Matmn(/(5)) (A = (aj^) G Matm„(I(5))) is a (not necessarily 
square) interval matrix, then the matrices A = (aj^) and A = (a^j) are called lower 
and upper matrices of the interval matrix A. 

Proposition 9 Let S he an idempotent semiring with a zero and a unity 1. The 
mapping A G Mat„„(/(5)) [A, A] G /(Mat„„(5)) is an isomorphism of idempotent 
semirings Mat„„(/(5)) anc? /(Mat„„ (5)). If the semiring S has an interval extension 
1(5), then this assertion remains true if I{S) is substituted by 1(5). 
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Here intervals [A, A] in /(Mat„n(S')) or I(Mat„n(5')) are defined with respect to 
the partial ordering ^ in Mat„„(5') (see example p.3|) . The proof follows easily from 
proposition ^ indeed, this proposition implies that addition and multiplication of 
interval matrices are reduced to separate addition and multiplication of their lower 
and upper matrices. 

The following proposition is an immediate consequence of theorem 0: 

Proposition 10 // a commutative idempotent semiring S with a zero and a unity 
1 is algebraically closed and satisfies cancellation and stabilization conditions, then 
for any matrix A G Mat„„(I(5')) there exist a nonzero 'eigenvector' \ G Matni(I(>S')) 
and an 'eigenvalue' [A, A] G 1(5*) such that AV = [A, A] V. // the matrix A is 
irreducible, then the 'eigenvalue' [A, A] is determined uniquely. 

It follows from proposition ^ that A V = A V and A V = A V. 
Consider the following interval discrete stationary Bellman equation (see also sec- 
tions ^, refs:matrices): 

X = AX0B, 

where A G Mat„„(I(S')), B,X G Mat„s(I(S')). Consider the following iterative pro- 



cess: 



Xfc+i = AXfc © B = A'^+iXo © A J B, (2) 

where X^ G Mat„,(I(S)), A; = 0, 1, . . . 

The following proposition is due to B. A. Carre [1^, Theorem 6.1] up to some 
terminology: 

Proposition 11 If matrix A G Matnn{S) is semi-definite, then the iterative process 
Xfc+i = AXk © B stabilizes to the solution X = A*B of the equation X = AX © B 
after at most n iterations for any initial approximation Xq G Mat„i(S') such that 
Xo 4 A*B. 

Suppose an idempotent semiring S satisfies the assumptions of proposition [1^. Let 
[Ai, Ai], . . . , [A^, Afe], 1 ^ k ^ n, he the eigenvalues of the matrix A G Mat„„(I(5')). 
Denote sup{Ai, . . . , A^} = A; by p(A). It is possible to give a simple spectral 

criterion of convergence of the iterative process (^: 

Theorem 4 Let S be a commutative semiring satisfying conditions of proposition [7^ 
and matrix A G Mat„„(I(S')) be such that p(A) ^ 1. Then the iterative process 
Xfc_|_i = AXfe © B, /c ^ 0, stabilizes to the minimal solution X = A*B of equation 
X = AX©B after at most n iterations for any Xq G Mat„i(I(S')) such that Xq =^ X. 



Proof. It follows from proposition ^ that it is sufficient to prove that sequences of 
lower and upper matrices of {X^} converge separately. To this end, we shall show that 



the matrices A and A are semi-definite; then the result will follow from proposition ITT 
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Since a^j =^ a^j for all i,j, we need only to prove that A is semi-definite. First we 

shall prove this if A is irreducible. Using the expression for a unique eigenvalue of an 
irreducible matrix A in terms of cycle invariants |]TB| 



J 
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i = l,...,n 



'■1112 







©a 



<f>{n)/l 



where (l){n) is the least common multiple of the numbers 1, . . . , n, we see that for any 
cycle its cycle invariant P = a^j^ © ■ ■ ■ a^^ satisfies P ^ 1 if A ^ 1 (indeed, if 

P © 1 ^ 1, then, by remark 0, (1 © P)'^(")/' = 1 © P'^(")/' ^ 1, so 1 © A"^^"^ y 1— a 
contradiction). Thus A is a semi-definite matrix. 

If A is reducible, there exists a permutation matrix Q such that A = QBQ^^, 
where B = (6ij)ij=i,...,n has an upper block triangular form: 



B 



/Pi * 
Pa 

\0 



* \ 



BkJ 



1 < k ^ n, 



and all square matrices Pi, . . . , P^ are either zero or irreducible. We claim that every 
eigenvalue A of A is an eigenvalue of some P;, I = 1, . . . ,k. Indeed, let V be an 
eigenvector of A with an eigenvalue A; denote ith element of the vector V by fj. 
Consider a decomposition {1, . . . ,n} = Ni U ■ ■ ■ U N^, where Ni (1 Nj = if i ^ j 
and P; = {bij)ij^Ni, I = I, ■ ■ ■ , k; let Iq = max{ / | 3i G Nf. Vi ^ 0}. We see that A 
is a unique eigenvalue of the irreducible matrix Bi^ corresponding to the eigenvector 
{'Vi)i(zNi^- The condition p(A) ^ 1 implies that Pi,...,Pfc are semi-definite. Since 
P = -< 1 for any cycle containing indices i E Ni, j E N^, I s, v/e conclude that A 
is a semi-definite matrix. □ 

Remark 12.1. Compare this simple proposition with the well-known spectral 
convergence criterion of the iterative process in traditional Interval Analysis ( |I2| , 
theorem 12.1]), which in our notation has the following form: 

The iterative process X^+i = AX^ + B, k > 0, converges to a unique solution X of 
the equation X = AX + B for any Xq G Mat„i(/(C)) if and only z/p(| A|) < 1. 
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